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ABSTRACT 

A new boundary condition is presented for simulating the flow over passively porous 
surfaces. The model builds on the prior work of R.H. Bush to eliminate the need for 
constructing grid within an underlying plenum, thereby simplifying the numerical modeling of 
passively porous flow control systems and reducing computation cost. Code experts for two 
structured-grid flow solvers , TLNS3D and CFUIX and one unstructured solver . USMTDns , 
collaborated with an experimental porosity expert to develop the model and implement it into 
their respective codes. Results presented for the three codes on a slender forebody with 
circumferential porosity and a wing with leading-edge porosity demonstrate a good agreement 
with experimental data and a remarkable ability to predict the aggregate aerodynamic effects 
of surface porosity with a simple boundary condition. 


INTRODUCTION 


Passive Porosity Technology (PassPorT) [ 1 ] 
is an enabling How alteration concept that can 
potentially resolve many aerodynamic problems. Its 
underlying principle is illustrated in the upper sketch ot 
Fig. 1 that depicts a porous skin positioned over a 
closed cavity/plenum region. Local pressure 
differences within flow over the outer surface 
“communicate” through the plenum in concert with 
small amounts of mass transfer in and out of the 
porous surface to alter its effective aerodynamic shape. 
For a properly designed system, the hole size is small 
with respect to the boundary layer thickness and is less 
than or equal to the skin thickness, and the flow 
velocity into and out of the plenum is low. PassPorT 
was originally applied to transonic airfoils to reduce 
the normal shock strength and thus, eliminate shock- 
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Figure 1 - The PassPorT concept and application 
conformal control effector on military aircraft. 




induced separation and lower the drag levels 1 2-7 ]. 

More recently the concept was applied to a 
supersonic delta wing to reduce the crossflow, shock- 
induced separation [8]. During this study, it was 
discovered that connecting the plenum to high and low r 
pressure regions on the wing, as conveyed in the 
middle sketch in Fig. I, allowed for localized “lift 
dumping” thus, providing roll and yaw control. 
Subsequently, many applications were envisioned such 
as conformal control effectors on military aircraft as 
illustrated on the lower half of Fig. I. The potential 
for broad application of this technology to 
aerodynamic concepts provided enough incentive to 
fund the fabrication of a number of porous models. 

A family of four Tangent-Ogive forebodies 
was built (two, 2.5 Caliber and two, 5.0-Caliber, each 
slenderness ratio consisted of one solid and one with 
22 r /( porosity surface). These models were tested in 
the NASA LaRC 7- by 10-Foot High Speed Wind 
Tunnel, 14- by 22-Foot Subsonic Wind Tunnel, and 
Unitary Plan Wind Tunnel. In addition, the 5.0- 
Caliber porous model was tested in the McDonnell 
Douglas Shear Flow 7 Facility. The results from these 
tests |9,10J showed that passive porous systems could 
be used to eliminate asymmetric loading conditions 
that occur at high angles of attack on axisymmelric 
slender bodies even at zero sideslip angle. This side 
force is a result of asymmetric pressure loading caused 
by differences in the location of the separation-induced 
vortex on either side of the forebody. This 
phenomenon has been documented extensively [II- 
1 4 1 . Numerous studies have been conducted to 
develop devices that eliminate or minimize the 
asymmetric behavior of forebodies [15-18). These 
"fixes" typically involve reshaping the nose or adding 
devices, such as strakes, to the existing geometry that 
add weight and are beneficial for only a limited range 
of conditions. In references 9 and 10, a passive porous 
system was investigated to alleviate the asymmetric 
loading on the forebody. These tests also showed that 
pitch and yaw control could be attained by opening or 
closing the porous surface on one side or the other on 
the forebodies. 

A zero-sweep, porous, GA(W)-1 [General 
Aviation (Whitcomb) - 1) [19] wing model was built 
and tested in the NASA LaRC 14- by 22-Foot 
Subsonic Wind Tunnel. All surfaces on the model were 
constructed with porosity that could be closed by 
covering the holes with tape. Through selective 
opening and closing of porous regions, the locations on 
the wing that induced the most pitch, yaw, and roll 
control were identified, as well as, the locations 
responsible for either increased performance (i.e., 
increased L/D nm and angle of attack that L/D llm occurs) 


or increased drag. This model was later reskinned and 
modified to include actuators inside the plenums to 
open and close the porous surfaces. The results showed 
that no pressure lag occurred. The recorded forces 
responded as quickly as the porous surface could be 
actuated indicating no appreciable lag due to pressure 
equalization. 

With the success achieved in the wind-tunnel 
tests, a series of research efforts were conducted to 
develop CFD models to represent a passive porous 
system. The first studies utilized Darcy's Law [20|. 
These had limited success and were fairly accurate as 
long as the coefficient required in the equation was 
chosen correctly. Later attempts with modified 
versions of Darcy's Law and utilizing some techniques 
used by researchers to determine oxygen transport 
through capillary walls were slightly more successful 
[21J. 

The work presented in this paper builds on a 
methodology first developed by McDonnell Douglas to 
simulate normal flow 7 through a screen positioned at a 
zonal-grid interface boundary [22). The primary 
contribution of this work is to reformulate that 
approach into a new 7 boundary condition for simulating 
the flow over porous aerodynamic surfaces. This 
eliminates the need for constructing grid within an 
underlying plenum, thereby simplifying the numerical 
modeling of passively porous flow control systems 
and reducing computation cost. Code experts for two 
structured-grid flow solvers. TLNS3D and CFL3D. 
and one unstructured solver. USM3Dns. collaborated 
with an experimental porosity expert to develop the 
model and implement it into their respective codes. 
This paper describes the formulation of the new 
boundary condition and presents an assessment of its 
effectiveness in simulating the aggregate aerodynamics 
induced by surface porosity using wind-tunnel results 
for a 5.0 Caliber Tangent-Ogive body and GA(W)-| 
wing model. This technique was recently utilized in a 
passive porosity control effector design study [ 1 1 on 
the military aircraft configuration illustrated in Fig. I. 

NOMENCLATURE 

A Area 

a Speed of sound 

b Wing span 

C n Drag force coefficient 

CFD Computational Fluid Dynamics 

C L Lift force coefficient 

C Pressure coefficient 

c Wing chord 

c, Magnification parameter 



c\ Relaxation parameter 

D Model base diameter 

UD Lift to Drag ratio 

M Mach number 

hi Area averaged mass flux 

l\f Number of cell faces on porous boundary 

p Pressure 

Re Reynolds number 

.v Solidity parameter, A svl JA j 

T Temperature 

US Upper surface 

it Local velocity normal to porous surface 

a\y,c Cartesian coordinates 

v Spanwise distance 

a Angle of attack (degrees) 

(p Contraction coefficient, or Circumferential 

station on slender body, 0 & 360 deg. - 
windward, 1 80 deg. - leeward 
y Ratio of specific heats for air, 1 .4 

p Density 

Subscripts 

avg averaged over surface 
h boundary 

c wing reference chord 

D model base reference diameter 

domain computational domain 
effective constricted orifice area 
max maximum 

min minimum 

normal perpendicular to porous surface 

/f iteration number 

porous region of surface porosity 
solid solid region of porous surface 

t total stagnation quantity 

tangent tangent to porous surface 
update iterative update quantity 

00 freestream condition 

1 upstream flow on porous surface 

2 minimum area of flow through porous surface 

3 downstream, fully mixed flow' through porous 
surface 

POROUS SURFACE FLOW MODEL 

The proposed porous surface flow model 
builds on the general approach of Bush [22]. The Bush 
model was derived to pass flow' information across a 
conterminous boundary separating an external 
computational grid zone from an interna! plenum grid 
zone. The present approach is formulated as a surface 


boundary condition, thus eliminating the need for grid 
cells on the plenum side of a porous surface. 

A schematic of the porous surface model is 
depicted in Fig. 2 as a cut through a section of solid 
surface and a hole. Porosity is defined by a solidity 
parameter s—A sifU /A i that quantifies the area ratio 
between solid and total surface. Three uniform states 
are assumed across the porous surface: upstream (I), 
minimum area (2), and downstream (fully mixed) (3). 
The upstream and downstream areas are equal {A =A,). 
The area at 2 is assumed to be reduced by the solidity 
parameter, s. and a contraction coefficient that 
represents a further area reduction due to flow through 
the orifice, A : =A,<p(l-s) where i^A^/A,). The 
contraction coefficient in Ref [22] is determined from 
curve fits to experimental data following guidelines 
provided by Cornell [23 1 an<i Rouse |24|: 

<p n = 0.04 1 37 /( 1 .0982 -(!->)) + 0.57293 + 0.005786( 1 - a > 

and 

V = <(>„+ 0.1 85 * 1 1/4 (/>,://>:-!)• 



Figure 2 - Crosscut schematic of porous surface. 


Governing equations 

The governing equations are derived from 
conservation of mass and momentum (normal to the 
surface) for steady, one dimensional, isentropic flow' 
of a perfect gas. Density and velocity are 
nondimensionali/ed by freestream density p ^ and 

speed of sound d„ , respectively, pressure by p^a^ « 
and temperature by static freesteam temperature 7^ . 

It is also assumed that the flow is adiabatic (i.e., no 
heat transfer), and that there is no total pressure loss 
during jet formation. The latter assumption is valid as 
long as the fluid passing through the orifice has an 
inviscid core. Thus, the flow normal to the porous 
surface boundary is determined by the following 
equations. 

Conservation of mass from 1-2: 

(/>!/), ={pll) 2 <p(\-.s) ( 1 ) 
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Conservation of mass from I -3 

(/?"), =(pu K 


( 2 ) 


No loss in jet formation from 1-2 
Pr\=Pt2 (3) 

Momentum balance from 2-3 with induced screen 
losses 

/':l' + yM\<p{\-s)\= />,(l + T/Ml) (4) 


From this result compute the maximum mass How 
using Eq. 5 written at state 3 


ipu )nu i \ 'Jr" ^ 


y - 1 > 

1 + ^ Mi 

i - 


( 8 ) 


If (pit ) , < (p«) m;lx . i. e. not choked, then solve Eq. 2 
using Eq. 5 written at stale 3 for M l : 


Auxiliary isentropic How relations: 

/ V/- 




/’, = p 
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pM 


I 


\ + l——M 2 

2 

\r/ (r-0 


(5) 


( 6 ) 


M 2 _ ~ 1 + a / 1 + 2 ((t- <)/r k {pu)j/ p* 

Y - 1 

then solve Eq. 1 using Eq. 5 written at stale 2. and Eq. 
4 for M 7 via Newton iteration: 


Assuming that a tangential pressure gradient 
exists across the surface and the net mass tlux across 
a passively porous surface is zero, there will always 
be some flow transversing both into and out of the 
plenum at any given time. Thus, the application of the 
governing equations is dependent on the direction of 
the surface normal How. which is a function of the 
local pressure gradient between the plenum and the 
domain. Formulations are presented below for the two 
conditions. 


Mi 


y-l O 

1 + * M ~ 

1 + yMj 

1 T 1V1 -> 

2 

1 + yM l<p{\ - s) 


<p 2 (\-sf 


t, (Hr 

r Py 


( 10 ) 


= 0 


If How is choked, (pu ), > (pu ) , then j 


set 


(pu)\ =(pu) mLn and Ml =1, Now use Eqs. 3 and 4 
with 6 written at slate 2 to compute p i{ : 


Implementation 

A general assumption for the present model is 
that entire process from state 1 to 3 in Fig. 2 transpires 
over an infinitesimal distance. Hence, the application 
of the model is in the form of a surface boundary 
condition. The following will describe the proposed 
boundary condition for flow moving both into {p vkmm < 

P, > and °u> of <p tl > J the plenum. The local 

velocity, u , is defined as the component normal to the 
surface. 


Pt\ = Py 


1 + 

I + }M ^(p{ 1 - L v) 


1 + 


y-1 


\r / (r- 1) 


M; 


(ID 


then compute M j 2 using Eq. 5 written at state 1 via. 
Newton iteration: 


Mr 


i + 


Y- 1 


M[ 


s 

/-I 


T, (pu)] 


= 0 


( 12 ) 


7' P7\ 

We now have sufficient information to compute the 
boundary flux using (/?//), , A/, 2 and p ir 


Flow r into plenum (p . < p , ) : 

* 1 J plenum r tltmuiw' 


The configuration for (low into the plenum 
from the domain is presented in Fig. 2. The first step 
is to assume choked How by assuming j and 

solve Eq. 1 using Eq. 5 written at state 3, and Eq. 4 
for M ; via Newton iteration 


Y + 1 


l + tfp(l-.s) 


<p 2 (\-s) 2 -M~\ 


+ 


Y - 1 
? 


Ml 



(7) 


Ph = P,\ 


i + ^1a/ 


r-i 


Pb = 


YPb 


1 + 


Y~ 1 


Mi 


\ 


) 


u h. normal = (p u ) I 1 Pb 

u h.hm*e*i zero * or a viscous boundary or the in viscid 
tangential velocity for an inviscid boundary. 
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Flow leaving plenum (/> > p iU „ m J 


The configuration for flow into the domain 
from the plenum is presented in Fig. 3. The first step 
is to once again check for choked flow by assuming 
- — i and solving Eq. 1 using Eq. 5 written at 

state 1, and Eq. 3 for M f via Newton iteration: 


r + 1 


i Y v- ] 

1+*^ — M; 


JJ 


(p 2 (\ -s) 2 -Mf =0 (13) 


From this result compute the maximum mass flow 
using Eq. 5 written at state 1 : 


. \ 


. 1/2 


Mf 


(14) 


If (pif), <(pu) nuy . i.e. not choked, then solve Eq. 2 
using Eq. 5 written ul state 1 lor M f : 




- ,+ tl 

i 

1 + 2' , 1, C 

V 

r Pi 


(15) 


r - 1 


then solve Eq. 1 using Eq. 5 written at stale 2. and Eq. 
3 lor M V via Newton iteration: 


1 + 


Y~ I 


M, 


2 

T, (p«)j 

r /t 




A7; 




r+i 

H) 


f a 0-*) 2 


= 0 


(16) 


If How is choked, (pit), > (p«) niax • then set 
(p (/ ) 3 =(pif) nux and Ml= 1. Now Eqs. 3 and 6 to 
compute pressure at state 2: 



Figure 3 - Model for flow leaving plenum. 


Pi = P\ 


( y - 1 > 

i + — — m r 






(17) 


then solve Eq. 2 using Eq. 5 written at stale 3 for M , . 


Ml 


1 + 


r-i 


Mi 


I + yM j<p( 1 - s ) 


1 + 


7*, _ 


(18) 


= 0 


r / ? 5 

Use Eq. 6 written at state 3 and Eq. 4 to compute p : 


Pit ~ Pi 


l + }M;(p{\ z s)( ] + Y Z l_ M 2 

1 ■ 


7 

y-\ 


(19) 


I + yMy 

There is now sufficient information to compute the 
boundary flux using (pit) } . and p f} . 


Ph = P,t 


1 + 


V 


Ph 


Wh 


1 + 


/-I 

2 

7 - 1 




A*i 


/-i , 


l hju>niuil ~ (P U )} f Ph 

tangential velocity for an in viscid boundary. 


u is zero for a viscous boundary or the inviscid 

h.uu itfi'irf 


Determination of Plenum Pressure 

A plenum pressure must be determined 
which yields a zero 1 net mass flux balance across the 
porous surface. This pressure is used tor p ? in Eqns. 

8-1 E and for p x in Eqns. 14-17. The following 
describes an iterative procedure for estimating 
plenum pressure. 

First compute an area-averaged pressure and 
normal mass flux over the porous surface. 


Pavfi 





a 


1 One could conceivably prescribe a non-zero net 
mass flux to simulate blowing or suction. This has not 
been tested in present work. 
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/h rnm.\ = 

Assume the initial value of plenum pressure to he the 
averaged pressure , = />„,,)• 

On the second iteration, compute a bracketed 
update to plenum pressure: 

P update ~ m ^ n ((l ( I pttrous )/ ? pU‘nm)Ln-\ ' Pt .° ■= ) 

/»^ f = ,nax (/'»,.</>»<■ ■' (H)5 />n,in ) 

where c ] =10 is a magnification parameter, p too the 

freestream total pressure, and p inm the minimum value 
of pressure on porous surface from summation for 

P a\ a ' 

It is necessary to impose a filter on the 
pressure updates to damp the temporal oscillations. 

P pfcnum.n ~ 2 P plenum. n-\ ~^Pupilatt‘ )/(* 2 0 

where the relaxation parameter r ■> = 50 . 


DESCRIP TI ON OF FLOW SOLVHRS 

The porous boundary condition has been 
implemented in two structured-grid flow solvers, 
TLNS3D and CFL3D, and one unstructured solver, 
USM3Dns. A brief description of the salient features 
of each code is included below. 

TLNS3D 

The TLNS3D (Thin Layer Navier-Stokes 
Solver for 3-D flows) code (25,26] uses a cell- 
centered, finite-volume approach for solving inviscid 
and viscous flows over complex configurations on 
general multi-block structured grids. A suite of 
algebraic and one- and two-equation turbulence 
models is used for simulating turbulent flows. 
Artificial dissipation is added to the central-difference 
scheme for stability. TLNS3D also makes use of local 
time stepping, implicit residual smoothing and 
multigrid techniques in conjunction with multi-stage 
Runge-Kutta time-stepping scheme to accelerate the 
convergence of the code to steady-state solutions. 

CFL3D 

CFL3D [27] solves the three-dimensional, 
time-dependent, thin-layer approximation to the 


Reynolds- Averaged Navier-Stokes (RANS) equations 
using a finite volume formulation in generalized 
coordinates. It uses upwind-biased spatial differencing 
with Roe's (28) flux-difference splitting or Van Leer’s 
[29| flux-vector splitting methods for the inviscid 
terms, and central differences for the viscous and heat 
transfer terms. The code, which is second-order 
accurate in space, is advanced in time with an implicit 
three-factor approximate factorization (AF) scheme. 
Temporal subiterations with multigrid are used to 
recover time accuracy lost as a result of the AF 
approach during unsteady calculations. The code 
includes several grid connection strategies, a vast array 
of zero-, one-, and two-equation turbulence models 
(linear as well as nonlinear), and numerous boundary 
conditions. The results presented in this paper were 
obtained using the two-equation k-to SST model of 
Menter |30|. 

USM3Dns 

USM3Dns [31] is a tetrahedral cell-centered, 
finite volume Euler and Navier-Stokes flow solver. 
Inviscid flux quantities are computed across each cell 
face using Roe’s 1 28] flux-difference splitting. Spatial 
discretization is accomplished by a novel 
reconstruction process that is based on an analytical 
formulation for computing solution gradients within 
tetrahedral cells. The solution is advanced to a steady 
state condition by an implicit backward-Euler time- 
stepping scheme. Flow turbulence effects are modeled 
by the Spalart-Allmaras one-equation turbulence 
model [321. The model can be integrated all the way to 
the wall, or can be coupled with a wall function to 
reduce the number of cells in the sublayer 

RESULTS AND DISCUSSION 

An assessment of the new r porous boundary 
condition is made with two selected examples using 
the three flow solvers. The first is a 5.0 caliber 
tangent-ogive with and without forebody porosity. 
Without porosity, the forebody develops strong vortex 
asymmetries at high angles of attack that generate 
severe yawing moments. Surface porosity alleviates 
these asymmetries. 

The second case is the General Aviation 
(Whitcomb) - I wing (GA(W)-l) with and without 
leading-edge porosity. The effect of the porosity is to 
reduce lift and increase drag. Such porosity is 
conceived as a conformal control device when applied 
asymmetrically to a wing to generate rolling and 
yawing moments. 
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5.0 Caliber Tangent-Ogive Forcbody 

A 5.0 caliber tangent-ogive forebody was 
tested, as one of a family of ogives, in the Langley 
7X10 Foot High-Speed Wind Tunnel to investigate the 
effect of fineness ratio on the asymmetric loading of 
slender forebodies and the etlectiveness oi passive 
porosity in alleviating these asymmetries 1 10]. Figure 4 
presents a full surface representation of the forebody. 
The original sting-mounted wind-tunnel model is 30 
inches long and 4 inches diameter. For the 
computational geometry depicted in Figure 4. the sting 
has been removed and the body extended to 40 inches 
in length where it is terminated at an outflow 
boundary. Surface porosity is applied to the darkened 
region w hich extends from a- I to 20 inches. 

The computation of asymmetric vortex Hows 
on slender ogive bodies at high angles of attack 
presents a challenging problem [33.34] that is beyond 
the scope of this paper. Such flows are characterized 
by massive crossflow separation w’ith asymmetric 
feedback through the boundary layer, and are highly 
sensitive to laminar-to-turbulent transition location, 
turbulence models, and numerical discretization. The 
present work is focused on the formulation and 
verification of the porous boundary condition. Since 
surface porosity restores flow symmetry, most of the 
following assessments will be performed on symmetric 
half-plane grids. One full-body asymmetric flow 
solution will be included ior completeness. 

Both structured-hexahedral and unstructured - 
tetrahedral grids were constructed for the study. The 
structured grid for TLNS3D and CFL3D contained 
177X97X65 cells, and the unstructured grid for 
USM3Dns had 1,009,929 ceils. The USM3Dns grid is 
mirrored for the full-body computation. Farfield 
boundaries were placed at x HnH , y wrw , and z mw of - 1 20, 0, 



Figure 4 - Surface representation of 5.0 
caliber tangent ogive configuration. Porous region 
denoted by shading. 


and -150 inches, and x Hwn . y„ tin , and z mn of 40, 180, 
and 150 inches, respectively. A characteristic 
inflow/outflow boundary condition was applied to the 
inflow', top and side boundaries. An outflow 
extrapolation condition was prescribed to the aft 
boundary. The no-slip condition was applied to the 
solid surfaces on the ogive, with the exception that 
USM3Dns utilized a wall function. A 22-percent 
porosity is =0.78) boundary condition is applied to 
the darkened forebody surface region denoted in Fig. 
4. 

Navier-Stokes solutions were computed on 
the ogive body with the three flow' codes at M 0.3. 
a=3() and 40 degrees, and Re { =()A million. TLNS3D 
and USM3Dns utilized the Spalart-Allmaras one- 
equation turbulence model and CFL3D the Menter 
SST two-equation model. The solid solutions were 
computed first. Then the porous boundary condition 
was activated and the solutions restarted. Figure 5 
presents the typical convergence ol the surface net 
mass flux, and the plenum and surface averaged 
pressures for the porous region after restarting from a 
converged solid surface solution. Note that the 
plenum pressure is substantially different than the 
surface-averaged value. Earlier studies assessed the 
use of surface-averaged pressure as the plenum 
pressure with the porosity boundary condition. This 
approach was determined to yield inaccurate results, 
which led to the development of the iterative 
procedure for plenum pressure described earlier. 



Iteration 

Figure 5 — Convergence of the net mass flux, and 
the plenum and surface averaged pressures over 
the porous region of 5.0-caliber tangent-ogive 
forebody from USM3Dns; M„=0.3, a=30 deg, and 
Re n =0.4 million. 
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Implementation of this procedure was only completed 
for USM3Dns in time for publication. For the results 
presented herein, the values for plenum pressure were 
determined from the USM3Dns solutions and 
provided as input for the other codes. The final 
nondimensional plenum pressures for a=30 and 40 
degrees were 0.71 18 and 0.7099. respectively, where 
nondimensional free-stream pressure is 0.7143. 
Nondimensionalization is discussed in earlier section 
on Governing Equations. 

The effect of porosity on the offbody How 
field is evident in Fig. 6. Contours of density are 
shown at station x-10 inches (a//7=2.5) for a=30 deg 
for the solid-surface solution (left) and porous- surface 
solution (right). The tightly clustered contour lines 
over the solid surface denote the presence of a strong 
vortex core, whereas the lines over the porous surface 
suggest a more diffused vortical system. 

A comparison of surface circumferential C 
distributions is shown in Figs. 7 and 8 at a— 10 inches 
( v//>=2.5) for angles of attack 30- and 40-degrees, 
respectively. A strong asymmetry is evident in the 
experimental distributions for the solid surface, which 
is the source of unwanted load asymmetries on the 
ogive. The beneficial effect of passive porosity is 
evidenced in the companion data by a diffusion of the 
vortex suction peaks into a symmetric distribution. 

A full-grid solution from USM3Dns is 
included in Figs. 7 and 8 to illustrate the difficulty of 
computing accurate surface pressure distributions with 
a solid (non-porous) surface for this class of problem. 
Note that asymmetric solutions are produced, but 
accuracy is marginal. These solutions were 



Figure 6. - Contours of density at x/D= 2.5. (Solid 
surface - left; porous surface - right.) 5.0 caliber 
tangent-ogive forebody from USM3Dns; M^=0.3, 
a=30 deg, and Re, =0.4 million. 


generated at zero sideslip with an initial asymmetry 
triggered by applying asymmetric viscous boundary 
conditions on the nose ahead of a=1; a wall function on 
the starboard and a no-slip condition on the port. 
Subsequent porous computations (not shown) were 



0 60 120 180 240 300 360 


<p, deg 

Figure 7 - Code to code comparison of 

circumferential C p distribution on 5.0 caliber 
tangent-ogive forebody with and without porosity. 
Station x/D=2.5, M^O.3, a=30 deg, and Re, =0.4 
million. 



Figure 8 - Code to code comparison of 

circumferential C p distribution on 5.0 caliber 
tangent ogive forebody with and without porosity. 
Station x/D=2.5, \L=0.3, a=40 deg, and Re, =0.4 
million. 
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restarted from the asymmetric solutions. Flow 
symmetry was restored even with the asymmetric 
viscous boundary conditions set ahead of ,v=l . 

The computational C, distributions on the 
half-plane grids are presented on Figs. 7 and 8 between 
<p = 0 and 180 degrees The distributions from the solid 
(non-porous) solutions are included to iacilitate code- 
to-eode comparisons, and are not intended to reflect a 
correct solution to an inherently asymmetric problem. 
Some differences in C /t distribution are noted between 
codes, but each qualitatively captures the dominant 
vortex How features. The locus is on the porous 
surface results which demonstrate good quantitative 
agreement with the experimental data. This confirms 
the adequacy of the new porous surface boundary 
condition for computing aggregate aerodynamic effects 
of passive porosity tor this class of problem. 

GA(W)-1 Wing 

The GA(W)-I [General Aviation (Whitcomb) 
- 1[ wing was tested in the 14- by 22-Fool Subsonic 
Wind Tunnel. This model had a 9-foot span, 3-fool 
chord, and 0° leading-edge sweep. Surface pressure 
taps were located at 3 spanwise locations. The model 
was also equipped with porous skins to represent a 
passive porous lest article (porous surface with plenum 
cavity). 

The semispan surface definition used for the 
computations is shown in Fig. 9. Porosity was applied 
to the shaded leading-edge region ahead of the 18- 
percenl chord station. Computational grids were 
constructed for each of the flow solvers. The grid lor 
TLNS3D an CFL3D contained 193X65X33 hexahedral 
cells, whereas for USM3Dns the grid contained 
1.681,831 tetrahedral cells. Farfield boundaries were 
placed 10 chord lengths away from the wing in all 
directions on which a characteristic inflow/outflow 
boundary condition was applied. The no-slip condition 
was applied to the wing solid surfaces, with the 
exception that USM3Dns utilized a wall function. 
When applying leading-edge porosity, a 22-percent 
condition (s~0.78) was prescribed to the darkened 
region denoted in Fig. 9. 

Navier- Stokes flow solutions were computed 
at A/«f 0.2, a=() and 8 degrees, and a chord Reynolds 
number of 3.5 million. As before, the plenum pressure 
was determined from USM3Dns and provided as input 
for the other codes. The nondimensional plenum 
values used for a=() and 8 degrees were 0.7130 and 
0.7156, respectively. 

Figures 10 compares the chord wise Cp 
distributions at a=8 deg. and 2y/b=0.67 (one chord 


length from the symmetry plane) between the code 
results and experimental data for the solid and porous 
surfaces. The experimental data reveals a dramatic loss 
of leading-edge suction peak and consequent loss of 
lift due to passive porosity. The solid surface 
computational results are nearly identical between the 
codes and are in generally good agreement with the 
experimental data, with the exception of the leading- 
edge suction peak. The porous leading-edge 
computations show' some variation between codes, but 
are in reasonably good agreement with the data. 

Figures 11 and 12 illustrate the large impact 
of leading-edge porosity on lilt and drag coefficients. 



Figure 9 - Semispan surface geometry for GA(W)-1 
wing. Porosity applied to shaded region around 
leading edge 



x/c 


Figure 10 - Cp comparison of solid and porous 
GA(W)-1 wing for TLNS3D and USM3Dns with 
experiment at 2y/b=0.67. Moo=0.17, a=8 deg., and 
Re =3.5 million. 
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and demonstrate that the porosity boundary condition 
model yields correct estimates of those effects at 
angles of attack of 0 and 8 degrees. With leading-edge 
porosity having such a large effect on lift and drag, an 
asymmetric application of this device could be 
envisioned for lateral-directional control in place of 
moving control surfaces. The computational model 
presented herein should be useful as a supplemental 
design tool in what was previously an experimental 
intensive process. 



-5 0 5 1 0 1 5 20 25 

a, deg 


Figure 11 - Comparison of lift characteristics of 
solid and porous leading edge for GA(W)-1 wing. 
M«=0.17 and Re ( = 3.5 million. 
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Figure 12 - Comparison of drag characteristics of 
solid and porous leading edge for GA(W)-1 wing. 
M w =0.17 and Re=3.5 million. 


CONCLUSIONS 

A new boundary condition is presented for 
simulating the flow over porous surfaces. The model 
builds on the prior work of Bush to eliminate the need 
for constructing grid within an underlying plenum, 
thereby simplifying the numerical modeling of 
passively porous How control systems and reducing 
computation cost. 

Code experts for two structured-grid How 
solvers, TLNS3D and CFL3D, and one unstructured 
solver, USM3Dns, collaborated with an experimental 
porosity expert to develop the model and implement it 
into their respective codes. Results presented for the 
three codes on a slender forebody with porosity, and a 
wing with leading-edge porosity demonstrate good 
agreement with experimental data and a remarkable 
ability to predict the aggregate aerodynamic effects of 
surface porosity with a simple boundary condition. 

Experimental studies of surface porosity have 
shown the strong potential for this technology as a 
flow control device. Porosity has many potential 
applications for aerodynamic control, drag 
reduction/production, separation control, and lift 
improvement. The present work should facilitate a 
more complete understanding of surface porosity in the 
future by enabling complementary computational 
studies and more timely design trades. 
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